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Abstract: Several work, is completed, within the 
framework of the modeling of a glow discharges 
continuous, one particular, the fluid, kinetic model and 
the hybrid model which is the combination between the 
simulation of Monte Carlo and the fluid model. In this 
paper, one uses the Monte Carlo method, to simulate a 
glow discharge continuous, in an electropositive gas 
(argon), has a pressure of 1Torr. This technique able to 
couple movement of particles charged with the 
Poisson's equation, it east makes it possible to calculate 
the spatial distribution of the electronic and ionic 
currents, the electric field distribution and also the 
Probability Density Function of electronic energies in 
plasma. 


Résumé : Plusieurs travaux, sont réalisé, dans le cadre 
de la modélisation des décharges luminescentes 
continues, en particulier, le modèle fluide, cinétique et 
le modèle hybride qui est la combinaison entre la 
simulation de Monte Carlo et le modèle fluide. Dans ce 
travail, on utilise seulement la méthode de Monte Carlo, 
pour simuler une décharge luminescente continue, dans 
un gaz électropositif (l'argon), a une pression de 1Torr. 
Cette technique capable de coupler le mouvement des 
particules chargées avec l'équation du Poisson, elle 
permet aussi de calculer la distribution spatiale des 
courants électroniques et ioniques, la distribution de 
champ électrique et également la fonction de densité de 
probabilité des énergies électroniques dans le plasma. 


Introduction: 


L’axe des travaux dans les décharges luminescentes 
s’est tourné vers une modélisation plus détaillée de la 
décharge, une bonne modélisation permet de 
comprendre et de prédire leur comportement général, 
alors de contrôler et d’utiliser cette technologie dans les 
besoins industriels, telque, les semi-conducteurs, les 
panneaux de visualisation à plasma, le dépôt des 
couches minces, la gravure en micro-électronique et le 
traitement de surfaces [1]. Il existe trois types de 
modèles, le modèle fluide qui consiste à résoudre les 
équations hydrodynamiques[2], [3], [4], ces dernières 
utilisent les moments de la fonction de distribution, qui 
sont accessibles à la mesure, mais ce modèle n’est plus 
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valable que dans le cas de l’hors équilibre 
thermodynamique, le modèle cinétique basé sur la 
résolution directe de l’équation de Boltzmann [5], [6], 
ou bien par un autre modèle statistique basé sur des lois 
stochastiques et simule les trajectoires des particules 
chargées,avec la technique appelée méthode de Monte 
Carlo[7],[8]. 


Description du modèle: 


Dans ce travail on utilise une géométrie cylindrique 
avec électrodes planes, le modèle à trois dimensions est 
réduit à deux dimensions axiale et radiale, dans ce 
modèle on considère qu’il n’y a pas d’effets radiaux 
(diffusion vers les électrodes), on considère aussi que le 
champ électrique est axial, et le modèle se réduit donc à 
une seule dimension, on garde la direction axiale et 
radiale de la vitesse. Le gaz de décharge est supposé à 
la température ambiante 300°k, on considère qu’il y’a 
quatre espèces des particules, les atomes neutres (Ar), 
les atomes excités (Ar”), les ions positifs (Ar*) et les 
électrons (e), la densité des particules chargées est très 
inférieure à celle des neutres. Pour cette raison on 
s’intéresse seulement aux collisions entre les particules 
chargées et les neutres, ces collisions sont: 

-Les collisions élastiques, d’excitation et d’ionisation 
entre les électrons et les atomes neutres, 
-Les collisions de transfert de charge et les collisions 
élastiques entre les ions et les atomes neutres de 
l’argon, les atomes excités se désexcitent spontanément. 


Simulation des particules chargées par la 
méthode de Monte Carlo: 

La méthode de Monte Carlo consiste à simuler un 
ensemble d’évènements (particules chargées) par un 
autre ensemble facilement réalisable (nombres 
aléatoires). Elle consiste à considérer une population 
des particules chargées et retracer l’histoire de chacune 
d’elle depuis son émission jusqu’à sa disparition, les 
trajectoires des particules chargées sont données par la 
loi de Newton. 


Z—Zo HvoAt+ 1 E.(AD? 


r=F +V,oAt (1) 


_ q 
Vz =Vot, EsoAt 
Vr —V,) 


zo et ro sont respectivement les positions axiale et 
radiale de la particule , z et r sont les positions axiale et 
radiale après un pas de temps Af; vs et v,osont 
respectivement les vitesses axiale et radiale à l’instant 
to. v.et v-les vitesses axiale et radiale à l’instant to+A, 
E est le champ électrique dans la position zo, et g et m 
sont respectivement la charge et la masse de la 
particule, l’énergie cinétique est donnée par la relation 
suivante: 


ba (2) 
2 r Zz 


Le champ électrique est obtenu par la résolution de 
l’équation de Poisson: 


2, 

= ie ge (ni) ne(2)) G) 
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ni(z) et ne(z) sont respectivement la densité des ions 
positifs et la densité des électrons. 

La probabilité de collision, la nouvelle énergie et la 
direction juste après la collision sont déterminées par un 
tirage des nombres aléatoires. 


P=I-exp(—-"m04(E)AS) (4) 
où P: la probabilité totale de collision, As est la distance 
entre le point (zoro) et le point (z,r), no la densité des 
neutres et (€) la section efficace totale de collision en 
fonction de l’énergie & .Le tirage d’un nombre aléatoire 
permet de dire si il y a une collision ou non, si R<P il y 
a une collision, si R>P pas de collision, maintenant si 
R<P, il faut tirer un autre nombre aléatoire (R ) pour 
déterminer le type de cette collision, en subdivisant 
l’intervalle des nombres aléatoires en sous intervalles 
chacun correspond à une fraction de probabilité, elles 
sont calculées à partir des relations suivants[9]: 

Pour les électrons: 


Oetoil €) =Oonis\ E)+Oexc(E)+O etat €) 


P er Oioni(£) 
ionis— 
Oero(E) 
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_ Oexc(E) 
BNC = ps 
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Pour les ions: 


(5) 


Oro E)=Oyan-charg(E) FO etasE) 


ee (CPR €) 


es arge (6) 
: hi Oiro(E) 


Oi,elas(E 
P: cas _ Sics(E) 


Oc.tot(E) 


OÙ Pionis Pexe Et Peelas SON respectivement les fraction de 
probabilité d’ionisation, d’excitation et des collisions 
élastiques avec les atomes neutres, Pan-charge €t Pielas 
sont respectivement les fractions de probabilité de 
transfert de charge et des collisions élastiques avec les 
atomes neutres. L’énergie de chaque particule et sa 
direction sont déterminées aussi par tirage des nombres 
aléatoires, après la collision la particule fait un parcours 
libre et la procédure se répète après un pas de temps. 


Simulation des électrons: 


Les électrons secondaires émis par la cathode sous 
l’effet du bombardement ionique, ces électrons ont des 
énergies distribuées entre 0 et 6 ev avec un maximum à 
4 ev, l’angle axial @ et radial g de chaque électron 
après l'émission sont données par les relations 
suivantes [9]. 
cosô, =R, 
po=2rR2 (7) 

Riet R2 sont des nombres aléatoires tirés au hasard, la 
nouvelle énergie et la direction dépendent du type de la 
collision. 

a)lonisation: dans le cas d’ionisation l’énergie total 
avant la collision est divisé entre l’électron primaire et 
l’électron secondaire. 


E0 — Eionis 
2 
E pri E0 —Eionis —Esec 


Esec =Ronis 
(8) 


OÙ &xoëmi et & sont respectivement l’énergie de 
l’électron secondaire, l’énergie de l’électron primaire et 
l’énergie de l’électron primaire juste avant la collision. 


Rionis: nombre aléatoire et &ons l’énergie d’ionisation 
d’argon  (£&onis=15.76 ev). 

b) Excitation: la nouvelle énergie est donnée par la 
relation suivante: 


E—Ep —Eexc (9) 
£exce et & Sont respectivement l’énergie d’excitation 
(&x—11.5 ev) et l’énergie de l’électron incident juste 
avant la collision. 
c) Collisions 
l’électron 

après une collision élastique est donnée par la relation 
suivante: 


élastiques: l'énergie cinétique de 


=a[1-2 Me ([_cos( ») (10) 


m, 


Où meet mx sont respectivement la masse de l’électron 
et la masse de l’atome cible et y l’angle de diffusion 
axial dans le repère de collision. 


fig. l:schéma représentant la séquence des collisions. 


L’angle axial juste après la collision est donné par la 
relation suivante: 
cos@=cos6, cos 7—sin 0, sin cosy (11) 


tel que ® l’angle axial juste avant la collision, 4 et w 
les angles de diffusion axial et azimutal. Ces angles 
sont calculés à partir de tirage des nombres aléatoires 
R;,et Ry. 


cos 7=I-R; 


W=27Ry (12) 
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* elastic collisions 
X ionization 
o excitation 


Sections efficaces en cm? 


Energie (ev) 


Fig.2. Section efficaces des collisions e’/ Ar [10], [11] 


Simulation des ions positifs: 


Dans le cas d’ionisation non seulement un électron est 
produit, mais aussi un ion positif, le transport de cet ion 
est simulé par la méthode de Monte Carlo, 
l’indépendance de simulation entre les électrons et les 
ions est due à la faible ionisation du gaz (les collisions 
électron-ion sont négligeables). L’énergie et la direction 
de l’ion juste après la collision dépendent du type de 
collision. 


a) Transfert de charge: dans ce type de collision la 
charge de l’ion se transfert vers l’atome et le nouveau 
ion démarre avec une vitesse nulle suivant la direction 
du champ, la section efficace de ce type de collision est 
supposée constante Oiran-charge =10°1* cm? [13], 


b) Collisions élastiques: les collisions entre les ions et 
les atomes d’argon sont considérées comme des 
collisions entre des sphères dures, la section efficace est 
donnée par la relation suivante: 
Oïelas =47Ra (13) 
Où R4r est le rayon de l’atome d’argon, l’énergie de 


l’ion juste après la collision est donnée par la relation 
suivante: 


Ejon=Ebiont 1-cos? 7) (14) 


OÙ &ion l'énergie de l’ion juste avant la collision et 7 
l’angle de diffusion axial dans le repère de collision. 


Résultats et discussion: 


A titre d'illustration sur l’utilisation de la méthode de 
Monte Carlo, on présente une application à une 
décharge continue luminescente pour un gaz d’argon à 
Itorr de pression, un espace inter-électrode de 2 cm et 
une tension continue sur l’électrode de gauche de (—300 
v) et l’autre est mise à la terre, le coefficient d’émission 
électronique secondaire est (y= 0.12). Les résultats 
numériques sont obtenus avec un pas temporel  Af = 
2.10"!isec et un pas spatial Az=0.04 cm. 

La figure (3) représente la densité des particules 
chargées (électrons et ion positifs), cette densité est très 
importante dans le plasma que dans les gaines car les 
électrons émis par la cathode sous l’effet de 
bombardement ionique gagnent l’énergie dans la gaine 
cathodique et excite le milieu d’une façon intense dans 
la lisière de la gaine. Le plasma est caractérisé par une 
neutralité électrique d’ou un champ électrique presque 
nul, par contre les gaines sont caractérisées par des 
charges d’espace positives qui sont responsables des 
chutes cathodique et anodique. 


+ Densité électronique 
* Densité ionique 


la densité en cm-3 


la distance en cm 


fig.3: la densité des particules chargées dans une 


décharge luminescente en continu 


Les figures (4) et (5) représentent respectivement, la 
variation spatial de champ électrique et la variation 
spatial de potentiel, la chute dans la gaine cathodique 
due à la densité des charges d’espace positives, dans le 
plasma, le champ électrique est pratiquement nul et le 
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potentiel est constant, car la densité des électrons est 
égale à celle des ions positifs, dans la gaine anodique la 
chute est due à la grande différence de mobilité entre 
les électrons et les ions et donc le potentiel du plasma 
est toujours supérieur à celle celui les électrodes. 


le champ électrique en (V/cem) 


In dintanan an nm 


fig.4: le champ électrique dans une 
décharge luminescente en continu 


le potentiel en (V) 


la distance en cm 


fig.5: le potentiel dans une décharge 
luminescente en continu 


Les figures (6) et (7) représentent respectivement la 
variation spatiale des densités des courants électronique 
et ionique et la variation spatiale de la densité de 
courant total, ce dernier est égal à la somme de deux 
courants ionique et électronique, dans la gaine anodique 
la densité des électrons est trop faible, de sorte que le 


courant total est presque complètement ionique, dans le 
plasma la densité des ions est importante, mais leur 
courant est faible à cause de leur faible mobilité et le 
courant total est un courant électronique, par contre 
dans la gaine anodique le courant électronique est 
supérieur à celui dans le plasma, ce dernier compense le 
courant ionique qui est du à la chute anodique, le 
courant total est pratiquement constant à tout endroit de 
la décharge avec une valeur de l’ordre de 0.0171 
mA/cm?. 


[ densité de courant ionique 
* densité de courant électronique 


la densité de courant en mA/cm2 


la distance en cm 


fig.6: le densité des courants électronique et ionique 
dans une décharge luminescente en continu 


la densité de courant total en mA/cm2 


la distance en cm 


fig.7: la densité de courant total dans une décharge 


luminescente en continu 
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La figure (8) représente la fonction densité de 
probabilité des énergies électroniques dans le plasma 
(FDP), elle est composée de deux parties avec une 
énergie moyenne de 6.29 eV , partie corps pour les 
faibles énergies où se trouve le grand nombre des 
électrons, ces derniers gouvernent pour l’essentiel les 
phénomènes de transports et la partie queue pour les 
énergies supérieures au premier seuil d’excitation, ils 
sont responsables des collisions inélastiques. 


FDP 


l'energie en (eV) 


fig.6: la fonction de densité de probabilité des énergies 
électroniques dans le plasma (FDP) 


Conclusion: 


L’élaboration d’un modèle statistique basé sur des lois 
stochastiques permet de simuler les trajectoires des 
particules chargées par une méthode de type Monte 
Carlo. Le couplage avec l’équation de Poisson donne 
un modèle auto-cohérent. De ce dernier, on peut prédire 
les caractéristiques générales de ce type de décharges et 
de déterminer d’une manière indirecte la fonction 
densité de probabilité des énergies électroniques, d’où 
la résolution de l’équation de Boltzmann, on peut aussi 
introduire facilement le second coefficient de Towsend, 
l'étude du comportement des électrons et des ions dans 
les situations hors équilibre thermodynamique (gaines). 
Alors, la prédiction et le contrôle du flux ionique et le 
potentiel au voisinage de la cathode, où les ions positifs 
sont accélérés et bombardent le substrat avec des 
énergies dépendent de la chute de potentiel dans la 
gaine, la distribution en énergie des ions détermine le 
dégât d’anisotropie de gravure. 
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